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FOREWORD 


The  Naval  Siirface  Warfare  Center  Dahlgren  Division  (NSWCDD)  has  a  growing 
responsibility  for  providing  ways  to  accurately  predict  and  detect  targets.  Information 
processed  through  track  filters  is  the  heart  of  all  major  combat  systems.  Tracking  and 
estimation  technologies  are  essential  to  the  future  of  the  United  States  Navy.  Currently, 
combat  system  functions  such  as,  Command  and  Decision,  Air  Control  and  Combat 
Systems  Component  Individual  Error  Checking,  rely  heavily  on  track  processing  data.  In 
the  future,  combat  system  functions  such  as.  Common  Command  and  Decision,  Multi- 
Platform  Weapons  Control,  Common  Air  Picture  and  Mission  Expansion,  will  rely  on 
track  processing  data. 

This  report  has  been  reviewed  at  NSWCDD  by  A.  Riedl,  Head,  Combat  Systems 
Technology  Division  (B30);  and  M.  Kuchinski,  Head,  Digital  Systems  Branch  (B32). 
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1  INTRODUCTION 


The  usage  of  traicking  filters  associated  with  a  tracking  radar  dates  back  to  work  by  Sklansky 
[20].  He  proposed  measures  of  performance  including  stability,  transient  response,  noise  and 
maneuver  error  as  a  function  of  the  dynamic  parameters  a  and  (3.  All  of  the  work  was  based  on  a 
frequency  domain  or  z-transform  analysis.  Subsequent  work  by  Benedict-Bordner  [2]  proposed 
a  relationship  between  a  and  ^  based  on  a  pole-matching  technique  that  combined  transient 
performance  and  noise  reduction  capability.  Subsequent  analysis  was  performed  by  Simpson 
[19] ,  Neal,  and  Benedict  [18]  for  the  a-p-j  filter.  By  this  time,  the  Kalman  filter  was  becoming 
well  known  in  the  radar  community.  Thereafter,  the  tendency  was  to  discuss  the  a  -  /3  and 
a- p-j  filters  as  steady  state  solution  to  the  Kalman  filter.  Subsequent  papers  by  Friedland 
[17],  Fitzgerald  [6],  and  Kalata  [13]  exploited  this  formal  similarity  to  derive  many  results  that 
can  be  used  to  characterize  tracking  performance  in  a  multi-tracking  environment.  The  basis 
for  the  analysis  of  performance  used  internally  with  the  Aegis  community  is  summarized  in  the 
internal  manual  entitled  “The  Working  Engineers  Guide  To  a  -  P  and  a  —  P  —  j  Filters”  by 
Reifler  and  Solomon  [23].  Later,  much  of  this  work  was  summarized  in  the  open  literature 
by  Kalata  [22].  A  summary  of  subsequent  developments  in  the  literature  to  1992  is  found  in 
Kalata  [22]  with  some  additional  work  since  then  found  in  Gray  ([10],  [11]),  as  well  as  the  open 
literature. 

Contrary  to  the  approach  that  is  usually  taken  in  the  literature,  we  propose  that  the  more 
‘natural’  viewpoint  is  to  introduce  the  constant  gain  filter  and  an  entity  that  is  independent 
of  the  Kalman  filter.  One  can  derive  the  information  that  characterized  the  filter  performance 
without  regard  to  the  Kalman  filter  design  criteria.  One  can  then  show  how  the  performance 
criteria  generalize  naturally  to  the  Kalman  filter  or  to  a  variation  of  the  Kalman  filter  that  re¬ 
places  the  process  noise  writh  a  bias  reduction  criteria.  While  variations  on  the  filters  discussed 
are  currently  in  use  within  naval  systems,  it  is  hkely  that  they  will  be  replaced  with  much  more 
advanced  estimation  techniques  such  as  the  interacting  multiple  model  (IMM)  filter.  Exploring 
the  ability  to  bound  filter  performance  is  a  necessary  part  of  the  redesign  to  replace  existing 
filters  with  advanced  filter  architectures.  An  IMM  design  that  consists  of  several  a-  P  and 
a-  p-j  filters  can  be  used  to  provide  such  a  bound  without  the  requirements  of  a  detailed 
system  simulation  [21] .  Different  architectures  can  be  explored  by  this  approach,  so  that  tight 
error  bounds  can  be  determined  as  part  of  an  overall  system  performance.  Subsequent  imple¬ 
mentation  of  a  true  IMM  would  then  be  used  and  known  to  have  performance  boundaries  within 
this  boundary.  We  will  take  the  results  derived  here  and  explore  such  issues  in  subsequent 
reports. 

The  a-  P  filter  has  found  application  when  large  numbers  of  objects  are  to  be  tracked.  By 
clever  selection  of  the  gains,  and  careful  design,  variable  gain  a  —  P  filters  combine  sufficient 
elements  of  the  Kalman  filter  ([8],[1])  so  that  there  is  not  significant  tracking  degradation.  Thus, 
there  is  useful  information  to  be  gained  by  a  detailed  performance  characterization  of  the  filter. 

The  tracking  equations  for  the  a  -  /?  filter  consist  of  two  parts:  prediction  equations,  which 
are  given  by 

Xp{k)  =  Xs{k  —  1)  +  Vs(k  —  i)T  (1-1) 
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Vp{k)  =Vs{k  —  1) 

and  smoothing  equations,  which  are  given  by 

(1-2) 

^s(k^  ~  Xp{k^  -|-  (x(Xffi(^k^  Xp{ky) 

(1-3) 

Vs(k)  =  vp{k)  -1-  ^{xm{k)  -  Xp{k)) 

(1-4) 

where 


•  Xs{k)  =  smoothed  position  at  the  k-th  interval 

•  Xp{k)  =  predicted  position  at  the  k-th  interval 

•  Xmik)  =  measured  position  at  the  k-th  interval 

•  Vs{k)  =  smoothed  velocity  at  the  k-th  interval 

•  Vp{k)  =  predicted  velocity  at  the  k-th  interval 

•  T  =  radar  update  interval  or  period 

•  a,  0  =  filter  weighing  coefficients 


Alternatively,  these  equations  can  be  written  as 


^13  +  G pXfji{l^ 


(1-5) 


where 


1-a  {l-a)T 


(1-6) 


ks)fc 


Xs{k) 

Vs{k) 


(1-7) 


(1-8) 


These  filter  equations  are  one-dimensional,  but  can  be  extended  to  three  dimensions  by  substi¬ 
tuting  successively  y  and  2  for  x  in  Eq.  (1-1)  through  Eq.  (1-4).  The  filter  equations  are  usually 
analyzed  in  one  dimension  and  the  resulting  analysis  is  usually  extended  to  three  dimensions 
with  the  assumption  that  similar  results  are  given. 


For  the  class  of  problems  when  this  occurs,  the  filter  can  be  viewed  as  a  constant  gain  filter 
which  is  nothing  more  than  a  matrix  difference  equation.  This  equation  can  then  be  solved 
regardless  of  the  measurement  model  provided  the  model  is  deterministic.  The  general  solution 
can  then  be  used  to  compute  the  covariance  matrix  under  very  general  assumptions  about  the 
noise.  This  is  an  alternative  and  slightly  more  general  to  work  done  by  Fitzgerald  [6], [7],  in  the 
early  eighties.  In  this  report,  the  general  case  will  be  solved  first,  and  the  filter  [2]  will  be 
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solved  as  an  illustrative  example.  Previously  [9],  the  z-transform  or  frequency  domain  method 
was  used,  but  here  the  direct  methods  that  have  become  more  fashionable  in  recent  years  will 
be  used.  The  solutions  are  independent  of  the  particular  relationship  between  a  and  P  that  are 
discussed  in  [2]  and  [13]. 

In  general,  the  update  equations  for  a  constant  gain  filter  can  be  written  as 

'Hn+l  =  +  Xn+lG,  (1-9) 


where 

vector  state  measurement, 

F,G  ^  filter  update,  gain  matrices, 

Xji  =  scalar  measurement  model. 

One  would  like  to  solve  Eq.  (1-9)  under  very  general  circumstances,  which  will  be  demonstrated 
in  the  next  section.  This  method  would  apply  to  any  scalar  measurement  model. 
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2  SOLUTIONS  TO  GENERALIZED  FILTER  EQUATIONS 


In  general,  the  update  equations  for  a  constant  gain  filter  can  be  written  as 


^n+l  =  FVn  +  Xn+lG,  (2-1) 

where, 

T]n  =  vector  state  measurement 
F,G  ^  filter  update,  gain  matrices, 

Xn  =  scalar  measurement  model. 

Hence,  for  n  =  0,  Eq.  (2-1)  has  the  form 


Vi  =  Fvo  + 

(2-2) 

For  n  =  1,  Eq.  (2-1)  has  the  form 

7/2  =  Frji  -f  X2G 

(2-3) 

Substituting  the  value  of  77^,  Eq.  (2-2),  into  Eq.  (2-3)  gives 

772  =  F{Fr)Q  +  XiG)  +  X2G 

(2-4) 

=  F^rjo  +  Fx\G  +  X2G 

For  n  =  2, 


%  =  -^’72  +  ^sG  (2-5) 

Substituting  the  value  of  772,  Eq.  (2-4),  into  Eq.  (2-5)  gives 

773  =  F{F^r}Q  -I-  FxiG  +  X2G)  +  X3G  (2-6) 

=  F^tjq F^xiG  +  FX2G -\- xsG 

and  so  on.  Thus  one  has  a  basis  for  induction,  and  can  then  establish  the  following  theorem: 

Theorem:  The  general  solution  to  Eq.  (2-1)  is 

k 

ri^  =  F'‘pQ  +  J2F''~"Gxn.  (2-7) 

n=l 

Proof:  Use  induction  or  substitution  to  verify  the  solution  directly,  as  done  above. 

One  can  recognize  that  in  the  theorem  the  general  solution  is  the  combination  of  the  homo¬ 
geneous  and  inhomogeneous  solutions.  The  homogeneous  solution,  found  for  Xn  equal  to  zero, 
is 


V'^  =  F%. 


(2-8) 
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The  inhomogeneous  solution  is 


k 


4  = 

72=1 


To  proceed  further,  the  powers  of  the  matrices  need  to  be  explicitly  evaluated. 
Let 

/(A)  =  det[F  -  AJ], 


by  the  Caley-Hamilton  theorem, 


/(F)  =  0  =  70-f  +  IxF  +  -  + 


(2-9) 


(2-10) 

(2-11) 


where  m  is  the  order  of  the  matrix.  Thus,  any  power  k  >  m  oi  a,  matrix  can  be  written  as 


F*'  =  7o(fc)^  +  7iik)F  +  ...  +  7m =  {F\  \j{k))  (2-12) 


where 


'  I  ■ 

7o(fc)  * 

F 

7i(fc) 

\F)  = 

7 

pm 

and  l7(fc))  = 

.  7m(*)  . 

(2-13) 


(note  {A\  =  \A)^  where  t  denotes  transpose.)  Note  the  solution  in  Eq.  (2-11)  can  be  simplified  by- 
expressing  as  in  Eq.(2-12)  provided  one  can  evaluate  sums  of  the  form  Yin  Therefore 

the  solution  is  of  the  form 


Vk  =  {F\  7(fc))  Vo  +  ^{F\  7(^  -  n))  Gxn- 


72=1 


(2-14) 
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3  a  —  (5  Filter 


The  a  —  /3  filter  has  found  application  when  large  numbers  of  objects  are  to  be  tracked. 
Thus,  there  is  useful  information  to  be  gained  by  a  detailed  performance  characterization  of  the 
filter.  The  tracking  equations  for  the  a  —  (3  filter  consists  of  two  parts:  prediction  equations, 
which  are  given  by 

Xp{k)  =  Xs{k  -  1)  +  Vs{k  -  1)T 
Vp(k)  =Vs{k-  1) 

and  smoothing  equations,  which  are  given  by 

““  ^p{k'^  “i“ 

Vs{k)  =  Vp{k)  +  ^{Xm{k)  -  Xp{k)) 

where 

•  Xs{k)  =  smoothed  position  at  the  k-th  interval 

•  Xp[k)  =  predicted  position  at  the  k-th  interval 

•  Xm{k)  =  measured  position  at  the  k-th  interval 

•  Vs{k)  =  smoothed  velocity  at  the  k-th  interval 

•  Vp{k)  =  predicted  velocity  at  the  k-th  interval 

•  T  =  radar  update  interval  or  period 

•  a,  0  =  filter  weighing  coefficients 

The  filter  gains,  a  and  0  satisfy  the  following  relation 

0  <  /?  <  a  <  1  (3-5) 

There  are  three  commonly  used  relationships  between  a  and  0.  The  first  is  the  Kalata  relation, 
which  is  obtained  from  steady  state  Kahnan  filter  theory  assuming  zero  mean  white  noise  in 
the  position  and  velocity  state  equations  [13]. 

0  =  2{2-a)-  (3-6) 

The  second  is  the  Benedict-Bordner  relation,  which  is  derived  based  on  good  noise  reduction 
and  good  tracking  through  maneuvers. 


(3-1) 

(3-2) 

(3-3) 

(3-4) 
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The  third  is  the  Continuous  White  Noise  (CTWN)  relation. 


Alternatively,  Eq.  (3-1)  through  Eq.  (3-4)  can  be  written  as 

l^s)fc  ~  l^p)fc  d"  GpXiji(k^ 


and 

where 


Xs{k) 

Vs{k) 


1_  \  _  I 

“  1  Up(A:) 


Qfi  — 


1  T 
0  1 

1-a  0 
1 

T  ^ 


Gp  = 


a 

§. 

T 


Alternatively,  Eq.  (3-9)  can  be  written  as 

l®s)fc  ~  ^j3  l^s)fc— 1  +  G  fiXm{k) 

where 

Fp  =  Hp-Qp: 

The  eigenvalues  of  Fp  can  be  shown  to  satisfy  the  equation 

/(A)  =  0  =  (1  -  a  -  A)(l  -  /3  -  A)  -t-  /3(1  -  a), 

which  simplifies  to 

/(A)  =  0  =  A^  +  (q^  +  /3  —  2) a  “h  (1  —  o:) 


1-a  (l‘-a)r 

-I  1-/? 


Let  r  =  yjl  —  a  and  2r  cos(^)  =  2  —  a  —  /3,  then 

/(A)  =  0  =  A^  —  2r  cos(0)A  + 


or 


/(A)  =  0  =  (A  -  re>^)iX  -  re"^’^). 


(3-8) 

(3-9) 

(3-10) 

(3-11) 

(3-12) 

(3-13) 

(3-14) 

(3-15) 

(3-16) 

(3-17) 

(3-18) 

(3-19) 

(3-20) 

(3-21) 
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Since  Fp  is  a  two  by  two  matrix,  it  can  be  written  as 

i^"  =  7oW+7i(fc)i"/3- 

The  eigenvalues  are 

A±  = 


(3-22) 


(3-23) 


Using  the  method  in  [5]  for  determining  the  power  of  a  matrix  gives  two  simultaneous  equations 
with  two  unknowns.  Define 

5(A)  =  7o(^)  +  7i(fc)A  (3-24) 


and 


When  fi(Ao)  =  ^(Ao),  one  gets 


h{X)  =  X’^. 


j-k^ke  _  ^ , 


and  when  fi(Ai)  =  g{Xi),  one  gets 


-je 


(3-25) 

(3-26) 

(3-27) 


Therefore,  there  exists  two  simultaneous  equations  to  solve  for  the  two  eigenvalues  to  get  the 
coefiicients.  Solving  these  two  equations  give  the  result 


7o(fc)  = 


— r*'sin((fe  — 1)6 


sin(5) 


and 


7i(fe)  = 


„fc-i 


sin(A:0) 


sin(5) 


(3-28) 


(3-29) 


The  general  solution  to  the  a  —  (3  filter  equations  consists  of  the  homogeneous  solutions  and 
inhomogeneous  equations.  As  stated  previously,  the  homogeneous  solution  is  of  the  form 


where 

therefore 


Vk=F%, 


F"  =  7o(fc)I  +  7i(fc)i"5> 


(3-30) 

(3-31) 


v1^  =  bo{k)I  +  'ri{k)Fi3]vo-  (3-32) 

By  substituting  the  known  values  for  7o(fe)  and  7i(A:),  the  homogeneous  solution  takes  the  form 


— 

■  Vk  = 


-r^  sin((fc  —  1)6)1  ^  sva.{k6)Ffi 

sin(5)  **”  sin(5) 


%• 


(3-33) 
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Also,  inhomogeneous  solutions  are  of  the  form 

k 

4  =  '£F''‘~"^nGp,  (3-34) 

n=l 


where 

pk-n  _  ^  (3-35) 

therefore 

k 

Vk  =  '^boik -n)I  +  ji{k-n)Fi3]xnGp,  (3-36) 

n=l 

fc-1 

=  J2boi^)I  +  ^li‘^)Ffibk-mGfi. 

m=0 

Again,  by  substituting  the  known  values  for  7o(fc)  and  7i(fc),  the  inhomogeneous  solution  takes 
the  form 

1  k—1 

4  =  ^  [— sin((m  -  1)6)1  +  sm{nW)Fp]xk-mGp.  (3-37) 

m=0 

One  can  further  simplify  the  inhomogeneous  solutions  in  Eq.  (3-37)  by  using  the  trigonometric 
identity 

sin[(m  -  1)0]  =  sin(m0)  cos(0)  -  cos(m0)  sin(0),  (3-38) 

therefore 


fc-i 


4=  E 


m=0 


/  (—  cos(0)  •  I  +  \ 

I  sin(m0)  -1-  (r”^  cos(m0))  1 1 


G/J. 


(3-39) 


Therefore  the  general  solution,  which  is  the  summation  of  homogeneous  and  inhomogeneous 
solutions  gives 

nk  =  [(-r''  sin((fc  -  1)0))  I  +  sin(A:0))  %  (3-40) 

fe-i 

+  E 

m=0 

where 

A  =  —  cos(0)7  +  — .  (3-41) 

r 


sin(0) 


sm{m6)  +  (r"^  cos 


(me))  I  Xk-n 


G/s, 
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3.1  CONSTANT  MOTION  MODEL 


If  one  assumes  that  Xk-m  is  a  constant,  say  xq,  then  the  general  solution  for  the  Constant 
Motion  Model  becomes 


nk  = 


sin(0) 


k—1 


sin((A:  —  1)0)  j  I  (r^  ^  sin(A;0)  j 


Vo 


(3-42) 


+  f  .  sin(m0)  +  (r’”  cos(m0))  I  ]  xqG/j. 

^0  Vsin(^)  J 


Examine  only  the  inhomogeneous  solution 


vi  =  sin(m0)  +  {r^  cos(m0))  Ij  xoGp 

xqGp 


(3-43) 


.  k-\  k-i 

■  ^  r’”sin(m0)  +  ^  r”^cos(m0)J 
^^^(  ^  m=0  m=0 


[sin(0) 


5(r,0)  +  G(r,0)/ 


XqG/S 


where 


and 


fe-i 

S{r,  0)  =  y]]  T’’"  sin(m0) 

m=0 

A;-l 

C{r,0)  =  ^  r^cos{m0). 

m=0 


Substituting  the  values  for  ^(r,  0)  and  C(r,  0)  (found  in  Appendix  A)  into  reveals 


^  ^  A  ( r  sin(0)  -1-  sin((fc  —  1)0)  —  sin(fc0) 

Vk  1  1  —  2r  cos(0) -I- 


sin(0) 
-h 


(3-44) 

(3-45) 

(3-46) 


1  —  r  cos(0)  -I-  cos{{k  —  1)0)  —  cos(fe0) 
1  —  2r  cos(0)  -1- 


/]xoG;3 


The  terms  in  the  inhomogeneous  solution,  Eq.  (3-46)  that  are  independent  of  k 


f  A  frsm{0)\  ,  fl-rcos{0)\  ^ 


simplify  to 


0 


(3-47) 


(3-48) 
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or 

xoU\,  (3-49) 

where 

«7i=  J  .  (3-50) 

The  inhomogeneous  solution,  Eq.  (3-46),  then  becomes 

rfk  =  xoUi  +  {^[-^^(r'^+^sm{{k-l)0)-r’‘sm{k9)^  (3-51) 

■f  cos((fc  —  1)0)  —  cos(A:0)^  I]xqGp}. 

Therefore,  by  substituting  the  simplified  inhomogeneous  solution  Eq.  (3-51)  into  Eq.  (3-42), 
the  general  solution  for  the  Constant  Motion  Model  is  of  the  form 

;i;^(^^sin((^- l)^K  +  ^''“^sin(A:0)E»?7o  (3-52) 

+UiXQ  +  sin((A:  -  1)0)  -  r’^  sin(A:0)) 

-I-  cos((fe  —  1)0)  -  cos(A:0)  j  IlxoGp}. 

Note  that  for  a  close  to  one,  the  filter  converges  quickly  to  steady  state.  While  for  smaller  a, 
the  filter  converges  very  slowly  to  steady  state. 
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Figure  3.1-3.  Constant  Motion  Model  with  a  =  .1 


3.2  LINEAR  MOTION  MODEL 


If  one  assumes  that  Xm  is  linear,  then  Xm  can  be  written  as,  Xk-m  =  k  —  m.  Thus,  the 
linear  is  separated  into  constant  and  linear  terms.  So  one  only  needs  to  calculate  the  linear 
term  and  then  apply  the  results  gained  from  the  constant  term  discussed  previously  in  Section 
3.1  to  get  the  complete  solution. 


^  sin(fc^))  E^]  t]q 


m 

■fc-i 


+ 


^0  (si^^""  sin(m0)  +  cos(m0)) 


m 


Examine  only  the  inhomogeneous  solution 

-A:-l 


^  sin(m^)  -t-  (r^  cos{m9))  A  m 

_m=0  ^  V  /  / 


G/q 


(3-53) 


(3-54) 
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k—l 


k-1 


mr‘^sm{m9)+  ^  (mr"^cos(m0))  J 

m=0  m=0 


Gp 


[sin(0) 


S,{r,e)  +  Ct{r,d)I 


Gi3 


where 


and 


fc-i 


Si{r,9)  =  mr'^sm{m9) 

m—O 


(3-55) 


k-1 


Ci{r,9)  =  y]]  mr'^cos{m9). 


(3-56) 


m=0 


Substituting  the  values  for  Si{r,9)  and  Ci(r,9)  (found  in  Appendix  A)  into  n|,,  reveals 


Vk  = 


A 


sin(0) 

A 


Siir, 9) +Ciir, 9)1^00 

-  ai  •  (a\  +  ((2  -  «)^  cos(6>)  -  2r^)  I)G^ 

sin{9)  ^  ^ 

+{  2^  /m  sin((fc  -  1)0)  -  (a  -  /3(1  -  A:))r*’  sin(fc0)) 

sin(0)  '' 

-h-i  ^(a  -h  kP)  cos((fc  -  1)0)  -  (a  -  p{l  -  k))  cos(fc0))}G^ 

r' 

The  terms  in  the  inhomogeneous  solution,  Eq.  (3-57),  that  are  independence  of  k 

A 


(3-57) 


/32  Lsin(0) 


ar  sin(0)  -I-  ^(2  —  a)r  cos(0)  —  I  G^, 


simplify  to 


a^/?  —  afS  -t-  a/?  —  ci^/? 
<Ai3  ,  ^ 

/ji  1  nn  /TT 


or 


Let 


T  T"  7^ 

0 
T 


i 

'T 


(3-58) 

(3-59) 

(3-60) 

(3-61) 


Therefore,  by  substituting  the  simplified  inhomogeneous  solution,  Eq.  (3-57)  and  the  non-k 
dependent  inhomogeneous  solution,  Eq.  (3-61)  into  the  general  solution,  Eq.  (3-53) 


i  _  T  ^  s  ^ 

~  0^^sm.{9) 


{^a  +  sin((A:  -  1)0)  -  (a  -  /3(1  -  fc))r*  sin(A:0)^  (3-62) 
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+  +  pk)  cos((fc  —  1)6)  —  {a  —  P{1  —  k))  cos{k6)'^}Gp. 

Combine  the  constant,  k,  with  the  linear  term. 

— m  “  ^  777-.  (3”G3) 

Therefore,  the  Linear  Motion  Model  solution  is  of  the  form 

Vk  =  ^j^^(^*'sin((fc-l)6>)/  +  r'^“^sin(fc6i)i=»77o  (3-64) 

+{kl  -  L) 

sin((A;  -  1)0)  -  (-a  +  /3)r*=  sin(fc0)) 

+  cos((fc  —  1)0)  —  (—a  +  /3)r*' cos(A:0)^  /]so<j/j}. 

Note  that  for  a  close  to  one,  the  filter  converges  quickly  to  steady  state.  While  for  smaller  a, 
the  filter  converges  very  slowly  to  steady  state. 
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Figure  3.2-3.  Linear  Motion  Model  with  a  =  .1 


3.3  QUADRATIC  MOTION  MODEL 


If  one  assumes  that  Xm  is  quadratic,  then  Xm  can  be  written  as  =  {k  —  rnf'  = 

—  2km  +  m^).  Thus,  the  quadratic  is  separated  into  constant,  linear,  and  quadratic  terms. 
So  one  only  needs  to  calculate  the  quadratic  term  and  then  apply  the  resiUts  gained  from  the 
constant  and  linear  terms  discussed  previously  to  get  the  complete  solution. 

^  sli^  ~  ^  sin(fe^))  F/j]  %  (3-65) 


+ 


sin(m^)  -I-  (r^  cos(m^)) 


m‘ 


Consider  only  the  inhomogeneous  solution  of  Eq.  (3-65), 

■fc-i 


Vk 


^0  ^  (r’”cos(m^))/^  m? 


(3-66) 
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•  4  fc-i  fc-i 

=  .  —  y~^  r'^  siTi{m9)m?  +  ^  {r'^  cos{mO))m^I 


fc-i 


A 


..  .  fc-i 

y~^  m^r'^  sm{m9)  +  ^  (rn?r'^  cos(Tn9)^  I 

m=0 


Gp 


sin(0) 


S2{t,9)+C2{t,9)I 


Gp 


where 


and 


fc-i 


S2{r,6)  =  ^  vn^r^  s\n{m6) 

m=0 


fc—l 


C2i'f’,9)  =  ^  Tn^r™^cos(7n0). 

m=0 

Substituting  the  value  of  S2ir,9)  and  C2(r,9)  (found  in  Appendix  A)  into  n|.,  reveals 
A 


(3-67) 


(3-68) 


4  = 


sin(0) 


S2ir,9)  +  C2{r,9)l\Gp 


(3-69) 


=  ^  +  a/3  —  2/3)r  sin(^)  —  Aar^I  -I-  {4a  —  2o?  —  a^')r  cos(0)  j')  Gp 

13^  Vsin(y)  / 

{(2a^  -20  +  al3  +  2aPk  -H  sin((A:  —  1)0) 


P^sm{9) 


-{2a‘^  -2/3  — a/3  +  2a/3k  + /3^  -  2k0^  +  k^0^)r^  sin(A:0)} 

-t-^{(2a^  -  2;0  4-  a/3  +  2a/3A:  -I-  k'^ j3‘^)r’^^^  cos((A;  -  1)0) 

-(2a^  -  2/3  -  a/3  -I-  2a^k  +  0^  -  2k0^  -h  k'^0^)r^  cos(A:0)}}G/3 
The  terms  of  the  inhomogeneous  solution,  Eq.  (3-69),  that  are  independent  of  k 


which  simplify  to 

0^ 

or 


4Q:r^/  +  (4a  —  2a^  -  a/3)r  cos(0)7j  Gp, 

(3-70) 

20{a  -  1)  ■ 

I3\2a  -  /3)  J  ’ 

(3-71) 

1 

-(2a-0) 

• 

(3-72) 

T/J 
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fci) 

^  —  -(2a-/3)  ’ 

T/3 

then  the  inhomogeneous  solution,  Eq.  (3-69),  becomes 

r]l  —  M  +  —  20  +  aP  +  2a/3k  +  k‘^0^)r'^'^^sm{{k  —  l)9) 

-{20?  -  20 -a(3  +  2a0k  +  0^  -  2k0^  +  k‘^0^)r^  sin(A:^)} 

+{{2a^  —  20  +  Oi0  +  2a0k  -I-  k'^ cos{{k  -  1)6) 

— (2a^  —  20  —  a0  +  2a0k  +  0^  —  2k0‘^  +  k‘^0^)0  cos{kd)}I}Gfi 


Combine  the  constant,  k"^,  and  linear,— 2/i:m,  terms  with  the  quadratic  term 

^k-m  =  —  2km  +  rr?. 

Therefore,  the  Quadratic  Motion  Model  solution  is  of  the  form 

nk  =  [(-r'=  sin((fc  -  1)6))  /  +  sin(fc6))  tjq 

+[{k‘^I  -  2kL  +  M) 

—  {2.0?  —  20  —  OL0  -f  ;0^)  sin(A:6)] 


-I-  ^2q;^  —  20  +  a/?)  cos{{k  —  1)6) J 
—  ^2q:^  —  20  —  a0  +  /3^)  r*'  cos{k6)I}G0 


(3-73) 


(3-74) 


(3-75) 


(3-76) 


The  terms  in  the  general  solution  of  the  Quadratic  Motion  Model,  Eq. 
transient  and  dependent  on  m 

-2kL  +  M 


(3-76)  that  are  not 
(3-77) 


simplify  to 


’  1  ‘ 
0 

-2k 

0 

1 

+ 

zlH-zS^.  1 

n  ( oc  l'\ 

\ 

\J 

V  . 

/ 

(3-78) 
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■  *="-2(V) 

2k-2(f-l) 

Recall 


so  one  can  define 


where  Ip  is  the  position  lag  in  the  response  due  to  the  acceleration  input,  and 


(3-79) 


(3-80) 

(3-81) 


4=  (2-^)  (M2) 

where  1^  is  the  velocity  lag  in  the  response  due  to  the  acceleration  input  [14].  The  lags 
mentioned  above  are  graphed  in  Figure  3.3-1. 


Since  0  <  r  <  1,  the  terms  and  represent  the  exponential  damping  in  the  transient 
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response.  Therefore,  the  non-transient  portion  of  the  smoothed  position  solution  is 

=  (3-83) 

and  the  non-transient  portion  of  the  smoothed  velocity  solution  is 

v,s  =  aoTk  -  2aolvT.  (3-84) 
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Figure  3.3-4.  Quadratic  Motion  Model  with  a  =  .1 


3.4  SINUSOIDAL  MOTION  MODEL 


If  one  assumes  that  Xm  is  sinusoidal,  then  Xm  can  be  written  as,  Xk-m  =  cos(A:  —  m),  then 
the  general  solution  for  the  Sinusoidal  Motion  Model  becomes 


^  sh^  ^  sin(A:0))  F>]  tjq 


(3-85) 


G, 


Note  that  for  a  close  to  one,  the  filter  converges  quickly  to  steady  state.  While  for  smaller  a, 
the  filter  converges  very  slowly  to  steady  state. 
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4  a  —  /?  —  7  FILTER 


The  tracking  equations  for  the  a  —  —  7  filter  consists  of  two  parts:  prediction  equations, 

which  are  given  by 

J.2 

(4-1) 


Xp{k)  =  Xs(k  -  1)  +  Vs{k  —  1)T  -f  -^as{k  -  1) 
Vp{k)  =  Vs{k  —  1)  Tas{k  —  1) 

Op(fc)  —as{k-  1) 

and  smoothing  equations,  which  are  given  by 


(4-2) 

(4-3) 


Xs{k)  =  Xp{k)  -1-  a{x.m{k)  -  Xp{k))  (4-4) 

Vs(k)  =  Vp{k)  -f  ^{xm{k)  -  Xp{k))  (4-5) 

as{k)  =  ap{k)  + -^{xmik)  -  Xp{k))  (4-6) 


where 

•  Xs{k)  =  smoothed  position  at  the  k-th  interval 

•  Xp{k)  =  predicted  position  at  the  k-th  interval 

•  Xra{k)  =  measured  position  at  the  k-th  interval 

•  Vs{k)  =  smoothed  velocity  at  the  k-th  interval 

•  Vp{k)  —  predicted  velocity  at  the  k-th  interval 

•  as{k)  =  smoothed  acceleration  at  the  k-th  interval 

•  GLp{k)  =  predicted  acceleration  at  the  k-th  interval 

•  T  =  radar  update  interval  or  period 

•  a, /3, 7  =  filter  weighing  coefficients 


For  the  a— /5— 7,  the  commonly  used  relationship  between  a  and  /3  is  the  Kalata  relationship 
which  is  obtained  from  steady  state  Kalman  filter  theory  assuming  zero  mean  white  noise  in 
the  position  and  velocity  state  equations  [13].  An  additional  relation  between  7  and  a  —  ^  is 
needed.  The  most  common  is 

0^ 


which  is  known  as  the  Neal-Simpson  relation[13]. 

Alternatively,  Eq.  (4-1)  through  Eq.  (4-6)  can  be  written  as 


(4-7) 


lxs)fc  =  \xp)f.  -H  GyXm{k) 


(4-8) 
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and 


where 


Q'i  i^s)/ 

k5)fc  = 

Xsik) 

Vs{k) 

as{k) 

\^p)k  ~ 

1 - 1 

1 _ J 

Q'i  — 


1 

0  1 
0  0 


rp  t2_ 
^  2 


T 

1 


1-a  0  0 


-f 


0 

1 


G/y  — 


a 

I 

T 

L  ^  J 


Alternatively,  Eq.  (4-8)  can  be  written  as 

\^s)k  ~  -^7  d" 


(4-9) 

(4-10) 

(4-11) 

(4-12) 

(4-13) 

(4-14) 

(4-15) 


where 

Fj  =  iJy  •  Q-y  = 

while  the  G  matrix  is 


1  —  a  (1  —  O')  T 


-T 

1-0 

—X 

T 

a 

G-y  = 

T 

.  ^  . 

If. 

2 


(1-0) 

(i-f)r 

^  2 


(4-16) 


(4-17) 


The  eigenvalues  of  Ky  can  be  shown  to  satisfy  the  equation  (using  the  result  in  Eq.  (4-16)  and 
expanding  about  the  33— minor) 

/(A)  =  0  =  (A  -  re^^){X  -  re-^^){l  -  |-A)  +  ^(3-;0-a)  +  (4-18) 

or 

/(A)  =  0  =  (A2  -  2r  cos(0)  +  r2)(l  -  |  -  A)  -  ^  (1  -h  2r  cos(0))  -  (^19) 
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which  simplifies  to 

/(A)  =  0  =  1  -  2rcos(0))  4-  A(r2  +  ^  cos(0))  - 

(4-20) 

/(A)  =  0  =  -A^b  +  A(b  +  r^  -s^'j  -r^ 

(4-21) 

where  7  =  1  -  and  -h  ^  -f-  2r  cos(0))  =  b.  Without  loss  of  generality,  one  can 

Eq.  (4-21)  can  be  written  as 

assume  that 

/(A)  =  0=(A-ie^'^)(A-fe-J>)(A-g) 

(4.22) 

=  0  =  (A^  —  2t A  cos{(p)  -1- 1^)  {A  — q) . 

Equating  Eq.  (4-21)  to  Eq.  (-^22)  leads  to  three  equations  in  three  unknowns: 

2t  cos  ip  +  q  =  b 

(4-23) 

-f  2tqcos(p  =  b  +  r‘^  — 

(4-24) 

qt^  = 

(4-25) 

Eq.  (4-23)  and  Eq.  (4-25)  have  solutions 

2tcos(p  =  b  —  q 

(4-26) 

J.2 

(4-27) 

while  the  Eq.  (4^24)  leads  to  the  equation 

—  (b  +  r^  —  -1-  br^'^  —  r^  =  0, 

(4-28) 

which  is  a  cubic  equation  in  t^.  Making  the  substitution  y  —  reduces  Eq.  (4-28) 

to  the  fundamental  form 

+py  +  m  =  0,  (4^29) 

where 

p=  --(b  +  r'^  -  -1- br^, 

(4-30) 

m  =  (b  +  r^-  ^  ~ 

(^31) 

Find  u  and  v  so  that 

3ui;  =  p, 

(^32) 

and 

—  v^  =  m. 

(4-33) 
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By  solving  for  u  in  Eq.  (4-32)  and  substituting  into  Eq.  (4-33),  the  result  is  given  as 


which  simplifies  to 


6  I  P  n 

V  +  mv  “  —  =  0. 


(4-34) 


(4-35) 


which  by  the  quadratic  formula,  has  the  solution 


3  =*=  V  ^ 


(4-36) 


To  continue,  one  needs  to  find  the  value  of  u  by  substituting  the  value  of  v  into  Eq.  (4-33) 


3  ~Tn  i  y  rn?‘  -|- 

«'='  =  m  H - - — 


(4-37) 


Hence,  the  value  of  y  can  be  found  by  subtracting  the  values  of  u  and  v  to  give  the  following; 


■  +  ^  3  —m  ±  Jm?  +  ^ 


(4.38) 


Therefore,  the  value  of  can  be  found  by  substituting  the  value  of  y  into  y  — 

^2  ^  -I-  ^  3  — m  ±  ^  ^  y.2  _  ^2 

2  \  2  + - 3 - 


(4-39) 


It  follows  that  the  values  of  q  and  q)  are  now  known. 

Let  h{X)  =  and  g{X)  =  7o(fc)  +  7i(^)'^  + 12{^)^-  By  equating  h(A)  and  5(A),  one  can 
see  that  if 

K^)  =  p(Ao)  (4.40) 


(4-41) 


^(-^l)  =  siM) 


(^42) 


and  if 


fk^  =  7Q(fc) -|- 7j(fc)te  -|-72(A:)t^e 


K^2)  =  9(^2) 


(4-43) 


(4-44) 


=  7o(^)  +  7i(^)9  +  J2ik)q^- 


(4-45) 
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Hence,  there  exists  three  simultaneous  equations  to  solve  for  the  three  eigenvalues  to  get  the 


coefl&cients[5] 

=  7q(A:)  -f  7i(fe)te^‘'^  -b 

(4-46) 

^  +  7i(fc)te-^V  +  72(fc)t2e-2^'^, 

(4-47) 

and 

=  7o(^)  +  7i(fc)9  +  72(fc)9^- 

(4-48) 

Solving  these  equations,  one  obtains 

7o(^)  =  C  sin  {{k  -  2)(p)  -  qH’^  sin  {{k  -  1)(,^))  -H  sin(¥?))  , 

(4-49) 

7^(fc)  =  ^  sin  ((fc  -  2)ip)  +  qt'^'’'-  sin  {(p)  -t-  qH  sin(2<,c>))  , 

(4-50) 

and 

72(A:)  =  ^  (t'^  sin  ((A;  -  !)</?)  -  qt’"~^  sin  {k(p)  +  /  sm{(p)^  , 

(4-51) 

where 

^  sin((/7)  [—2qt  cos{(p)  +<p'  + 

(^52) 

The  general  solution  to  the  a  —  /3  —  7  filter  equations  consists  of  the  homogeneous  solution 

(4-53) 

where  F'°  = 

jQ{k)I  +  'ri{k)Fy  +  ^2ik)F^,  hence 

Vk  =  +  7i(fe)^7  +  12{^)F^  Vo- 

(4r54) 

By  substituting  Eq.  (4-49)  through  Eq.  (4r52)  into  Eq.  (4-54),  the  homogeneous  solution  can 
be  written  as  follows: 

Vk  =  C{^gt^'*‘’'sin((fc  -  2)(/5)  —  g^t*^  sin((A:  —  1)</?) -l-9^t^sin(^)^  / 

(4-55) 

-h  sin  {{k  -  2)(p)  +  qt'^~^  sin  ((^)  -t-  qHsin{2(p)^  F^ 

+  (t'^  sin  {{k  —  l)(p)  -  qt^~^  sin  {k(p)  +  ^  sin(¥?)^  -F’7  }%- 

Recall,  the  inhomogeneous  solution  is  of  the  form 

rfk=j^F^-^XnG^,  (4-56) 

n=l 


where 


'^  =  7o(fc  -  n)I  +  —  n)F^  +  72(fc  -  n)F^ 


(4-57) 
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hence, 

k 

Vl  =  X)[7o(^-’^)-f  +  7i(^-^)-p7  +  72(^-”)-P’7]^n<?7  (4-58) 

n=l 

k—1 

=  S  [7o("^)-f  +  7l("i)i^7  +  72(l'n)F^]Xk-mG^. 

m=0 

By  substituting  Eq.  (4-49)  through  Eq.  (4r52)  into  Eq.  (4-58),  the  inhomogeneous  solutions 
can  be  written  as  follows; 


k-l 


Vk  =  ^  X)  ^  sin  ((m  -  2)^)  -  sin  ((m  -  1)^)  -j-  sin(v7) j  I 

(4-59) 

+  sin  ((m  -  2)(p)  +  qf^~^  sin  {(p)  -f  q^tsm{2(p)^ 

+  (f^ sin  ((m  -  1)9?)  -  qt‘^~'^  sin  (mp)  +  q^  sin(<^)^  F^}xk-mG^. 

One  can  further  simplify  the  inhomogeneous  solutions  in  Eq.  (4-59)  by  using  the  trigonometric 
identities 

sin[(m  —  2)p]  =  (cos^[</?]  —  sin^[^])  sin[m(/?]  —  (2sin[9p]  cos[^])  cos[my?] 

(4-60) 

and 

sin[(m  —  1)(^]  =  sin(77iy?)  cos(^)  —  cos{mp)  sin(9?). 

(4-61) 

therefore 

fc — 1 

ri  =  52  [((St’"sin(m6)  +  (Dt”^  cos(m0))  -f  E)  Gp, 

m=0 

(4.62) 

where 

B  =  -g(gcos[¥5]  - 1  cos[2y7])/  - 1  cos[2(,o]Ey  -1-  t{-q  -f  t  .cos[9?])E^, 

(4-63) 

D  =  sin[(/j]{q(q  -  2tcos[p])1 4-  2tcos[^]Ey  -  F^}, 

(4-64) 

and 

E  =  sin[v?]{9”^t^7  +  -|-  2qtcos[p])Fy  +  q^F^}. 

(4-65) 

Therefore,  the  combination  of  homogeneous,  Eq.  (4-55),  and  inhomogeneous,  Eq. 
solutions  lends  the  general  solution  in  the  following  form: 

(4-62), 

^ik  =  (9^^"'’^  sin  {{k  -  2)p)  -  qH'^  sin  ((A:  -  1)(^)  4-  qH^  sin(^)  j  I 

(4-66) 

4-  sin  {{k  -  2)p)  +  qt'^-^  sin  {p)  4-  sin(2^))  Fy 

4-  (t’^  sin  {{k  —  l)p)  —  qt^~^  sin  {kp)  4-  q^  sin(^)^  ^^7)^0 
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k-l 

+  ^  [{{Bt™'  sin(m^)  +  {Df^  cos(m<^))  +  E)  Xk-m]  Gp. 

m=0 


4.1  CONSTANT  ACCELERATION  MOTION  MODEL 


If  one  assumes  that  Xm  =  xq  +  VQm  +  ^oom?,  then  the  general  solution  for  the  Constant 
Acceleration  Motion  Model  is  of  the  form 

Uk  =  sin  ({k  -2)(p)-  qH^  sin  {{k  -  1)(^)  +  sin((/7)^  I  (4-67) 

-1-  sin  {{k  —  2)ip)  +  qt^~^  sin  (v?)  -I-  q^t  sin(2(^)^  Ey 

-f  sin  ((A:  -  1)^)  —  qt'‘~^  sin  (kcp)  -f  g*  sin(</3)  j  -^7  }% 
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The  figures  listed  below  of  the  Constant  Motion  Model  illiistrate  the  convergence  properties 
of  the  filter.  Note  that  for  a  close  to  one,  the  filter  converges  quickly  to  steady  state.  While 
for  smaller  a,  the  filter  converges  very  slowly  to  steady  state. 


Sfcnootii 

0- 

Predicted 

Function 

Time(Sec) 

Figure  4.1-1.  Constant  Acceleration  Motion  Model  with  a 


-  -0-  •  Predicted 
-  Rancticn 


Time(Sec) 

Figure  4.1-2.  Constant  Acceleration  Motion  Model  with  a 
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Figure  4.1-3.  Constant  Acceleration  Motion  Model  with  a  =  .1 


4.2  JERK  MOTION  MODEL 


If  one  assumes  that  Xm  =  ,  then  the  general  solution  for  the  Jerk  Motion  Model 


becomes 


Uk  =  sin  ((fe  -  2)<p)  -  qH^  sin  ((fe  -  l)(;o)  +  qh'^  sin(v?))  I 

sin  ((A:  -  2)ip)  -I-  sin  (v?)  -I-  sin(2^)J 
sin  ((fc  -  l)(p)  -  sin  {k(p)  -h  g*"  sin(<^))  ■F’7  }^o 


-k 

-t- 


fc-i  r 

+  E 

771—0 


sm.{Tmp)  -I-  {Df^  cos{m(p))  +  E)  ^  ^ 


Gp. 


Consider  only  the  inhomogeneous  solution  of  Eq.  (4-68), 


k-l 


rj\,  =  ^  ^  {{Bf^  sin(mv?)  -I-  {Df^  cos{rrnp))  -t-  E) 


/  yom^\ 


Gp 


(4-68) 


(4-69) 
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where 


and 


^yoT^ 


k—1 


E  f^sm(m(p)  +  cos{mip)^  +Em^^  Gp, 


m=0 


fc-i 


{BSzit,^)  +  {DCz{t,  <p))  +Y,Em^]  Gp 

m=0  / 

k-1 

Szlt^tp)  =  ^  sin{m(p), 
m=0 

fe--l 

Cz{t,(p)  =  mh'^  cos(m(p), 

m=0 

B  =  -q{q  cos[v?]  - 1  cos[2v3]) J  - 1  cos[2(^]F^  +  t{—q  + 1  cos[v3])F^ 

D  =  sin[(^]{5(5  -  2tcos[(p])I  +  2f  cos[<^]Fy  —  F?}, 


(^70) 

(4-71) 

(4-72) 
(4-73) 

E  =  -|-  -|-  2qtcos[(:p])F^  +  q^F^}.  (4-74) 

By  evaluating  S3{t,(p)  and  C3{t,ip),  as  done  in  Appendix  A,  and  considering  only  the  terms 
independent  of  m.  The  terms  independent  of  m  in  the  inhomogeneous  solution  simplify  to 


Recall 


which  implies 


where  Ip  is  the  position  lag, 


where  is  the  velocity  lag,  and 


where  la  is  the  acceleration  lag. 


(i-a)r^ 


Xs 

Vs  , 
as 


,  _  (l-a)r3 

jp 

^  r\/ 


(4-75) 


(4-76) 


(4-77) 


(4-78) 

(4-79) 
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1  2  3  4  5  6  7 

Time(Sec) 


Figure  4.2-3.  Jerk  Motion  Model  with  a  =  .1 


4.3  SINUSOIDAL  MOTION  MODEL 

If  one  assumes  that  Xm  =  sin(mT),  then  the  general  solution  for  the  Sinusoidal  Motion 
Model  becomes 

TTfc  =  sin  ((fc  —  2)</?)  —  g^t^sin((A;  —  1)(^)  +  sin(^)^  J  (4-80) 

-I-  sin  {{k  —  2)(^)  d-  sin  (c^)  -I-  c^t  sin(2^)  j  Ey 

-t-  sin  ((fe  -  1)^)  -  sin  (k(p)  -h  q’^  sin((/j)^  F^}Vo 

fc-i 

+  ^  sin{m(p)  +  {Df^  cos{mip))  -|-  E)  (sin(mr))]  Gp. 

m=0 
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The  figures  listed  below  of  the  Sinusoidal  Motion  Model  illustrate  the  convergence  properties 
of  the  filter.  Note  that  for  a  close  to  one,  the  filter  converges  quickly  to  steady  state.  While 
for  smaller  a,  the  filter  converges  very  slowly  to  steady  state. 


Time(Sec) 

Figure  4.3-1.  Sinusoidal  Motion  Model  with  a  =  .9 


Time(Sec) 

Figure  4.3-2.  Sinusoidal  Motion  Model  with  a:  =  .5 
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5  NOISE  REDUCTION  RATIOS 


Recall  that  one  can  express  the  update  equations  for  a  constant  gain  filter  in  the  following 
form 

Vk+i  =  +  ^k+iG,  (5-1) 

where  E{xn)  =  0  and  independent  of  t?*..  The  covariance  matrix  of  a  first  order  system  is 
defined  as  [15]  (where  t  denotes  transpose) 

Pk  =  E{{vk-EiVk)){Vk-E{Vk)n  (5-2) 

=  E{{nk  -  nk){nk  -  Uk)*, 

where  E(r)k)  =f}k  =  F’nk-v  One  can  calculate  Pk  by  substituting  Eq.  (5-1)  into  Eq.  (5-2) 

Pk  =  E{  {Frjk^i  -t-  XkG  -  Ffjk-i)  (-F’^l/t-i  -I-  ^kG  -  Ffik-i)*}  (5-3) 

=  E{Fr]k-iVk-iF*  +  4GG^}. 

The  expected  value  is  explicitly  given  as 

Pk  =  f2{P^k-iVk-iF'  +  4+iGG^}  (5-4) 

A:=0 

=  (E  ri,-A-,]  F'+if:  4+.)  GG* 

\fc=0  /  \fc=o  / 

=  FPk-iF^+<rlGG* 

where  E{xnXn)  =  4-  steady  state  {Pk  =  Pk-i  =  P),  so  the  covariance  is  in  the  form  of  a 
Lyapunov  matrix  equation 

Fppi  -p  =  -alGG\  (5-5)  . 

Define  a  matrix  S  such  that 

PE‘  -FP  =  S  (5-6) 

By  adding  Eq.  (5-5)  and  Eq.  (5-6),  one  can  solve  for  P  explicitly  in  terms  of  known  quantities 
to  give 

FPF*  -P  +  PF^  -FP  =  S-  0-2  GG*  (5-7) 

(F  +  I)P{F*  -I)  =  S-  alGG^  (5-8) 

P={F  +  I)-\S  -  o-2GG‘)(F‘  -  /)-!  (5-9) 

We  now  have  an  expression  of  P.  To  solve  for  S,  we  pre-multiply  Eq.  (5-6)  by  F  and  post- 
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multiply  Eq.  (5-6)  by  which  gives 

FP  (F*y  -  (F)2  PF^  =  FSF^  (5-10) 

then  substituting  from  Eq.  (5-5) 

FPF^  =  P  -  (tIGG^  (5-11) 

into  Eq.  (5-10)  gives  the  solution  for  S 

P5P*  -  5  =  FalGG^  -  alGG^F*  (5-12) 

One  can  then  solve  this  equation  for  S  and  then  substitute  S  into  Eq.  (5-9)  to  get  the  value  of 
P. 


5.1  NOISE  REDUCTION  RATIOS  FOR  THE  a- p  FILTER 
For  the  a  —  P  filter,  recall 


F3  = 


1-a  (l-a)r 


and 


a 

A 

T 


By  substituting  Fp  and  Gp  into  Eq.  (5-12) 

FpSF^^  -S=  [FpGpG^p  -  GpG^pF^^  <r: 

FpSFl  -S  =  a: 

Solving  for  S,  yields 


0  0^ 
-0^  0 


S  = 


//? 


0  p^ 
-0^  0 


where 


//3  =  -ap 

One  can  now  solve  for  P  by  substituting  S  into  Eq.  (5-9) 
P  =  {F^  +  I)-\S-<TlGpG^p){F^p-I)-^ 


(5-13) 

(5-14) 

(5-15) 

(5-16) 

(5-17) 

(5-18) 

(5-19) 


2-0  -T(1-q)  1 

4—20'—/?  4-'2a— /3 

/5  2-a 

■ 

-|(a  -  f) 

-1 

-r(l-a) 

1  1 
T 
—a 

r(4-2a-j0)  4-2a-0  J 

T 

L  p 

P  J 
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-2o?+0(Za-2)  -/3(2a-)3) 


a(-4+2Q+/3)  ar(-4+2a+/3) 

-0(2a-l3)  -20^ 

L  QT{-4+2a+/?)  ar-i(-4+2a+/3) 

The  elements  of  this  covariance  matrix  reveal  the  noise  reduction  ratios  for  position  (P®), 
velocity  {Pv),  and  position  x  velocity (Prt,). 


p  cn'l  -  ~ 

a(-4  +  2Q  +  /3)  ’ 


P;(0) 


-2/32 


ar2(-4  +  2a  +  /?)  ’ 


and 


-P{2a-P)_ 

-  ^j^(_4  +  2a  +  ^)' 


(5-20) 

(5-21) 

(5-22) 


5.2  NOISE  REDUCTION  RATIOS  FOR  THE  a -(3--/  FILTER 
For  the  a  —  /3  -  7  filter,  recall 


F^  = 


1-a  (l-a)T  (1-a)^" 


4 


l-(3 

T 


(i-f) 

**-  2 


2 

T 


and 


a 

T 


By  substituting  and  into  Eq.  (5-12) 

F-ySF^  -  5  =  (F^G^Gi,  -  G^G*^F*)  al 

n  2jg^-2ajg+/37  7(2/3+7) 

U  oT*  27^ 


-5  =  4 


-20^+2a0-0y 

2T 

— K(2i9+7) 

2r-^ 


2T 

0 


r 

0 


Solving  for  S,  yields 


"4 


0 

1 

1 

0 

7 

0  ■ 

0 

7^/5  +  ^ 

-h 

-7 

0 

0 

-7^5 

0 

0 

0 

0 

(5-23) 


(5-24) 

(5-25) 

(5-26) 


(5-27) 
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and 


Pa{0)  = 


_ _ _ 

T^(2al3  +  a7  -  27)(4  -  2a  - 13)  ’ 


P  rm-  2/37(2/3-7)4 

'  T^{2ap  +  a'y-2j){4-2a-p)' 


(5-39) 

(5-40) 
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6  COST  FUNCTIONS  AS  ALTERNATIVE  TO  TRACKING  INDEX 


Independent  of  which  of  the  three  relationships  between  a  and  P  one  assumes,  each  rela¬ 
tionship  can  be  shown  to  obey  the  common  constraint  due  to  Kalata  [13] 

r2  =  (6-1) 

1  -  a 

where  the  variable  T  is  commonly  known  in  the  naval  community  as  the  Kalata  tracking  index, 

r2  =  (6-2) 

^  m 

The  tracking  index  is  a  function  of  the  assumed  target  maneuverability  variance  al  (deviation 
from  modeled  behavior),  radar  measurement  noise  variance  and  T  is  the  update  interval. 
The  maneuverability  is  an  unknown  parameter  in  most  cases  because  there  is  no  direct  means 
of  determining  it  from  S3rstem  parameters,  nor  is  it  measurable.  One  invents  a  process  noise 
model  that  is  considered  the  best  means  of  modeling  unknown  threat  behavior.  The  variance 
in  maneuverability  is  cr„. 

There  are  alternative  means  of  selecting  the  values  for  the  filter  coefficients  than  the  tracking 
index.  This  method  is  based  on  the  fact  that  any  unmodeled  error  in  a  filter  introduces  a  bias 
which  is  characteristic  of  the  filter  response.  A  cost  function  that  combines  noise  reduction  and 

bias  is  introduced  and  minimized  with  respect  to  the  a  if  used  in  the  filter  equations  (one  has, 

of  course,  reduced  p  to  p  =  p  (a)).  This  methodology  amounts  to  accepting  a  mean  squared 
error  as  the  arbitrator  of  performance  characterization.  One  could  adapt  other  criteria  as  well 
[24],  but  mean  square  is  widely  accepted. 

The  response  of  an  a  -  filter  to  a  linear  acceleration  produces  three  terms;  a  transient 
term,  a  lag,  and  the  model  input  term.  The  lag  is  a  bias,  so  the  expected  value  E  of  the  steady 
state  response  E[xs{k)  -  Xm{k)]  is  the  lag.  Similarly, 

£;[(x.(fc)  -  xm{k))^]  =  Px{0)<tI  +  \LlalT^  (6-3) 

and 

E[{vs{k)  -  Xm{k)f]  =  PviO)ol  +  LlalT'^.  (6-4) 

This  allows  one  to  define  a  specific  cost  function  [23]  for  the  smoothed  velocity 

Jv{a,P)  =  E[{vs{k)  -  Xm{k))^]  (6-5) 

=  ajF.(0)  +  (2-i)"o?r^ 

where  the  terms  have  been  previously  defined.  By  introducing  a  change  of  variable 
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one  can  cast  the  second  term  of  the  cost  function  that  is  invariant  with  respect  to  the  relationship 
chosen  between  a  and  (3.  The  velocity  noise  reduction  ratio  can  be  written  as 


Pv{0) 


-2/?2 

ar2(-4  +  2a  +  ^) 


|T2(-4  +  2a  +  /?) 
rr2(4  -  20t  -  0) 

which  can  be  cast  purely  in  terms  of  r  once  the  relationship  between  a  and  0  is  made.  Therefore, 
the  specific  cost  function  for  smooth  velocity  is  of  the  form 


Jv{oc,0)  =  cr2p„(0)  + 


rT2(4  -  20r  -  0) 


For  the  Benedict-Bordner  relation  0qq  =  substituted  into  "t  =  §  gives 


(r  +  1) 


When  one  uses  the  Benedict-Bordner  filter  relation,  one  finds  the  velocity  noise  reduction  ratio 


= 

■*  V 


tT^{A-20t  -  0) 


(6-11) 


rT2(2r2-l) 

For  the  Kalata  relation  0  =  2(2  —  a)  —  4\/l  —  ot,  substituted  into  t  =  §  gives 


(2r  +  l)^ 


(6-12) 


0^  = 


{2r  +  lf 


(6-13) 
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When  one  uses  the  Kalata  filter  relation,  one  finds  the  velocity  noise  reduction  ratio 


yK 


2/3 

tT2(4  -  2/3r  - 13) 
4 

rr2(4r2  - 1) 


(6-14) 


For  the  CTWN  relation  a  =  ^^2/3  4- -  f ,  substituted  into  t  =  |  gives 

12t 


and 


/3^  = 


(6t2  4-  6t  -I-  1) 
12 


(6r2  4-  6r  4- 1) ' 

When  one  uses  the  CTWN  filter  relation,  one  finds  the  velocity  noise  reduction  ratio 

pC  _  2/3 

-  tT2(4-2/3t-/3) 


(6-15) 

(6-16) 

(6-17) 


tT^  (3r2  - 1) 


In  general,  when  one  uses  the  different  filter  relationships,  one  finds  the  velocity  noise 
reduction  ratio  (T=l) 

(6-18) 


■  ^  T  (<ir2  —  1) 

where  d  =  2  for  Benedict-Bordner,  d  =  4  for  Kalata,  and  d  =  3,  for  CWTN.  Eq.  (6-8)  can  now 


be  written  as  a  normalized  cost  function  (e.g.  =  J) 

^  n 


(6-19) 


where  Tr  =  .  Note  that  F/j  has  been  rewritten  with  a  subscript  to  avoid  confusion  with  the 

Kalata  tracking  index  which  has  the  same  physical  units  but  a  different  interpretation.  There  are 
a  couple  of  interesting  things  to  note  about  the  cost  function.  First,  the  bias  has  been  reduced 
to  a  form  that  is  invariant  with  respect  to  the  relationship  between  the  filter  coefficients.  This 
means  that  one  can  minimize  the  cost  function  to  arrive  at  r  =  T(r),  hence  one  can  arrive  at 
a  selection  criteria  that  has  an  absolute  invariance  with  respect  to  lag.  Additionally,  one  can 
maintain  the  measure  of  performance;  namely  that  the  filter  coefficient  are  minimized  jointly 
between  velocity  noise  reduction  and  “velocity  lag”.  This  leads  to  different  performance  than 
the  Kalata  criteria  for  selecting  the  filter  coefficients  which  is  based  on  a  plant  noise  model  as 
a  means  of  selecting  the  filter  coefficients. 
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To  compare  the  performance 
how  to  compute  the  coefficients, 
setting  it  equal  to  zero  gives 


of  different  coefficient  selection 
Taking  the  derivative  of  Eq. 


,2  _  d  (Sdr^  —  1) 

^  T^{dT^  —  1)2  (2r  -  1)  ’ 


techniques,  one  can  first  show 
(6-19)  with  respect  to  r  and 

(6-20) 


This  can  be  solved  numerically  or  graphically  to  give  r  =  r  (Tr)  which  solves  the  filter  selection 
coefficient  problem.  Substituting  Eq.  (6-20)  into  Eq.  (6-19)  gives 


J  optimized) 


d  d{3dT^  -  1)  (t  -  5) 

r  (dr2  -  1)  ^  2r2  {dr^  -  1)^ 


(6-21) 


which  is  cost  function  evaluated  at  the  optimized  value  of  r.  The  tracking  index  is 


T^  = 


r2  (t2  -  1) 


(6-22) 


for  the  Benedict-Bordner  relationship 


r2 


64 

16r‘^  —  8r2  —  l 


(6-23) 


for  the  Kalata  relationship, 


r2 


144 

Shr"*  —  24r2  -4- 1 


(6-24) 


for  the  CTWN  relationship.  To  make  a  comparison,  one  would  substitute  each  one  of  these 
equations  into  Eq.  (6-19)  and  compare  the  plot  for  all  r. 
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For  a  given  r,  one  sees  that  the  Kalata  relation  gives  the  smallest  a  for  a  given  r . 


Figure  6.1-5.  Isolated  Comparison  Between  a  and  r 


The  result  of  this  comparison  is  that  one  will  have  smaller  cost  for  the  optimization  technique, 
which  is  indicative  of  better  performance.  Note,  this  performance  improvement  occurs  when 
there  is  immodeled  behavior  only,  so  this  result  is  useful  in  a  complex  target  environment,  but 
not  necessarily  in  a  benign  threat  environment. 

Different  cost  functions  can  be  used  for  problems  other  than  the  ones  already  discussed 
might  prove  useful  under  some  circumstances.  If  there  is  concern  about  transient  effects  that 
occur  while  the  target  is  accelerating,  replacing  the  lags  by  twice  their  value  includes  this  effect 
since  the  transient  response  is  always  less  than  or  equal  to  the  lag.  Another  possibility  is 
to  weight  the  lag  by  the  percentage  of  time  maneuvers  are  expected  to  occur.  In  some  cases, 
it  is  more  desirable  to  minimize  the  lag  in  for  the  predicted  filter  response  rather  than  the 
acceleration  lag.  If  minimum  mean  square  in  position  is  the  primary  design  criteria  (say  as 
might  be  used  in  track  maintenance),  then  applying  the  same  analysis  for  the  position  cost 
function  is  appropriate.  Furthermore,  exactly  the  same  type  of  cost  function  approach  can  be 
applied  to  the  a  -  /?  —  7  filter.  Further  discussion  of  the  details  of  this  will  be  presented  in 
a  future  report  which  will  concentrate  on  combining  these  results  with  a  track  maintenance 
implementation  of  a  multiple  model  a  —  /?/q:  —  /?  —  7  filter  that  will  be  used  to  illustrate 
how  to  bound  IMM  performance  in  a  manner  useful  to  Naval  applications.  The  next  section 
illustrates  how  the  cost-function  technique  can  be  used  in  an  implementation  of  a  table-driven 
filter  implementation. 
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6.1  DESIGN  EXAMPLES 


To  illustrate  that,  with  proper  initialization  of  the  radar  track  data,  the  a  —  {3  filter  can 
be  implemented  as  a  look-up  table  of  the  coefficients  as  a  function  of  range  achieves  similar 
performance  as  a  Kalman  filter.  Unless  one  is  tracking  only  a  few  objects,  the  a  —  (3  filter 
implementation  of  the  Kalman  filter  could  be  preferred  over  the  conventional  Kalman  filter  for 
real-time  AAW  applications.  It  does  require  some  care  in  the  design  in  terms  of  xmderstanding 
the  operating  environment  as  well  as  careful  thought  in  the  underlying  design  concepts.  With 
that  understanding,  one  can  consider  several  specific  examples.  When  tracking  is  done  in 
Cartesian  coordinates,  the  noise  is  range  dependent.  Typically,  the  noise  can  be  written  as 
ax  =  Rag,  where  ag  is  the  sensor  angular  noise,  which  is  a  known  parameter  of  the  tracking 
system.  Given  that  the  range  is  in  kilometers  and  the  angle  noise  is  in  milliradians,  the  noise 
can  be  written  as  a  function  of  the  range  (n)  times  a  constant  k.  The  tracking  index 

kn  ■ 

Specific  system  parameters  are  then  plugged  into  Eq.  (6-25)  and  the  ex's  are  computed.  For 
example,  commercial  aircrafts  do  not  have  maneuvers  that  exceed  2g's.  For  simpHcity,  choose 
other  system  parameters  so  that  the  tracking  index  is  F  =  -1^.  Table  6.1.1  shows  the  two 
different  a's  computed  from  the  common  tracking  index  for  ranges  from  4  —  128  km. 


64 

.03 

.21 

.16 

128 

.014 

.15 

.12 

Table  6.1-1. 


Cost  Function  versus  Range 


OiR 

Pi^(O) 

.71 

.46 

.12 

.38 

.61 

.51 

.23 

.52 

.38 

.04 

.10 

.43 

.28 

.01 

.04 

.35 

.20 

.006 

.019 

.28 

.15 

.002 

.008 

Two  other  examples  of  interest  to  the  tracking  community  can  be  mapped  into  the  same 
form  as  the  accelerating  target  with  different  interpretations  of  the  coeflSicients  used  to  form  the 
tracking  index.  A  maneuvering  target  that  has  a  turning  rate  of  less  than  15°  deg/sec  has  an 
acceleration  that  is  equivalent  to  a  constant  acceleration  so  the  constant  in  the  tracking  index 
can  be  represented  by  where  cj  is  the  turning  rate  and  po  is  the  radius  of  ciuvattue  of  the 
turn.  The  noise  is  now  angular  an  ^  Rag  and  to  a  98%  confidence  level,  a  becomes  ZagR/T. 
Thus,  this  example  is  completely  mapped  into  the  solution  in  Eq.  (6-25)  with 


ZRag 


(6-26) 


Another  example  is  to  determine  the  ballistic  coefficient  of  an  object.  The  force  on  an  object 
undergoing  a  ballistic  slowdown  is  F  =  ^pkiP' .  This  force  produces  an  equivalent  acceleration 
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term  oq  = 


(-)" 


which  can  be  used  to  give  the  tracking  index 


,  pk'^T^ 

^  Zrri^VQRaQ ' 


(6-27) 


These  examples  illustrate  the  usefulness  of  this  approach  in  that  there  is  a  great  deal  of 
flexibility  to  attack  different  problems,  that  can  all  be  mapped  into  the  same  tracking  index 
with  a  different  definitions  of  its  parameters. 
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7  DETERMINATION  OF  FILTER  COEFFICIENTS 


The  selection  of  the  relationship  between  a  and  /?  remains  somewhat  problematic  in  the 
discussion  of  filters.  Because  of  the  close  connection  between  the  steady  state  solution  of  the 
Kalman  filter  and  the  constant  gain  a  -  /3  filter,  there  has  been  no  independent  development 
of  the  coefficient  relations  other  than  the  Benedict-Bordner  relationship.  This  can  be  rectified 
as  show  below,  but  it  remains  largely  an  academic  exercise  for  many.  Thus  it  may  be  skipped 
by  those  who  are  uninterested  without  serious  loss  of  continuity  with  the  previous  portions  of 
the  document. 

Recall  that  cost  functions  in  the  previous  chapter  are  used  to  determine  a  in  terms  of  known 
system  and  design  parameters.  This  cost  function  approach  can  also  be  used  to  determine  filter 
coefficient  relationships  ^  =  P(a).  Recall  from  optimization  theory  [?],  that  for  a  function  of 
two  variables,  say  x  and  y,  it  is  possible  to  minimize  this  function  so  that  a  minimum  exists  with 
respect  to  some  objective  criteria.  Once  the  relevant  equations  are  solved,  a  single  functional 
relationship  y  =  y{x)  is  achieved  that  satisfies  the  minimization  criteria.  To  illustrate  how  to 
do  this,  note  that  a  cost  function  J  is  minimized  by  taking  the  gradient  and  taking  the  inner 
product  with  a  vector  normal  to  that  function 

VJ  •  n  =  0.  (7-1) 


To  illustrate  this,  consider  the  fortuitous  choice  of  a  cost  function 

Taking  the  directional  derivative  oiJ{x,y)  with  respect  to  x  and  y  respectively,  gives  the  matrix 
equation  (with  the  normal  vector  corresponding  to  the  matrix  of  parameters  a  and  b) 


(fxh-fh^)  .  1 

r  -1 

h:^{x,y) 

a 

Uy^-fhy).  Q 

b 

h^{x^y)  yy 

=  0. 


(7-3) 


Clearly  this  equation  is  satisfied  only  if  the  determinant  is  zero,  which  gives 


9yifxh  ~  fhx)  —  dxifyh  fhy)  —  0- 


(7-4) 


Solving  this  then  implies  y  =  y{x),  which  is  equivalent  to  a  filter  coefficient  relationship. 

There  are  four  obvious  choices  for  the  cost  function.  We  have  already  discussed  two.  Recall 
that  the  noise  reduction  ratio  for  the  smoothed  velocity  is 


PviO) 


20^ 

T2a(4  -2a-p) 

W 

T2r(4  -  2r/?  -  0) 


(7-5) 
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so  a  velocity-lag-noise-reduction  cost  function  is 

Recall  that  the  noise  reduction  ratio  for  the  smoothed  position  is  given  by 

P(n)  =  2a^+/?(2-3a) 

^  a(4  -2a -13) 

2/3r^  +  (2  -  3tI3) 
t(4  —  2r/3  —  P) 

SO  a  position-lag-noise-reduction  cost  function  is 

2/3r2  +  (2-3T/3)  Al-r/S)^ 

JnLp  =  a — - - +  0- 


(7-6) 


r(4  —  2tI3  —  j3) 


(7-7) 


(7-8) 


Two  other  choices  for  cost  functions  date  back  to  Benedict  [2].  The  transient  velocity  perfor¬ 
mance  of  an  a  —  /?  is 


1  a2(2  -  a)  +  2y0(l  -  a) 

^  al3{A  -2a-p) 

1  I3t^{2  -  I3t)  -I-  2(1  -  (3t) 
r2  ^r(4  -  2r/?  -  0) 

SO  a  transient- velocity-noise-reduction  cost  function  is 


(7-9) 


Jntv  =  a- 


+  b 


/3t^(2  —  /3r)  +  2(1  —  /3t) 


(7-10) 


(7-11) 


t(4  —  2tI3  —  /?)'“  .  /3r(4  —  2t/3  —  /3) 

The  transient  position  performance  of  an  a  —  is 

^  (2  -  a)(l  - 

^  a;0(4  —  2a  —  0) 

(2-/?r)(l-^r)2 
t/?2(4  -  2^r  -  P) 

so  a  transient-position-noise-reduction  cost  function  is  (which  is  the  criteria  used  by  Benedict 
to  derive  the  BB  relationship,  though  he  used  a  somewhat  convoluted  z-transform  technique) 


_  2/5t2  +  (2  -  3r0)  (2  -  /?r)(l  -  /3r)2 

T{4-2rP-P)  '^\p\A-2Pr-py 


(7-12) 
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7.1  VELOCITY-LAG-NOISE^REDUCTION  COST  FUNCTION 


To  determine  which  different  coefficient  relationship  are  achieved,  first  consider  the  velocity- 


lag-noise^reduction  cost  function,  Jnlv-  Note  (from  now  on  there  is  the  mapping  x  = 

y  = 

f{r,P)  =  l3, 

T  and 

(7-13) 

which  implies 

/r(T,/3)  =  0, 

(7-14) 

and 

fl3{r,l3)  =  1. 

(7-15) 

Also, 

-  V  O 

II 

1 

to! 

b 

(7-16) 

which  implies 

Pr(r,/3)  =  (2t- 1), 

(7-17) 

and 

9p{r,  13)^0. 

(7-18) 

Finally,  note 

h(r,  /?)  =  t(4  -  2r/3  -  /?); 

(7-19) 

which  implies 

hT{r,p)  =  (4-4r/5-/?), 

(7-20) 

and 

=  -(2r-f  l)r 

(7-21) 

Additionally, 

(7-22) 

and 

-  fhp)  =  At 

(7-23) 

are  needed.  Substituting  Eq. 

(7-17),  Eq.  (7-18),  Eq.  (7-22),  and  Eq.  (7-23)  into  Eq. 

(7-4), 

gives 

4r  -  8t2  =  0 

which  implies  that  there  is  no  optimal  relationship  between  the  coefficients  for  this  cost  function. 
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7.2  POSITION-LAG-NOISE-REDUCTION  COST  FUNCTION 

To  determine  which  diSerent  coefficient  relationship  are  achieved,  first  consider  the  position- 
lag-noise-reduction  cost  function,  Jnlp',  note  that 


=  2/?r2  +  {2-  ZrfS), 

(7-24) 

which  implies 

=  4/3r  -  3/5, 

(7-25) 

and 

Also, 

ffi{T,P)  =  {2T^  -Zt). 

(7-26) 

5(t,/3)  =  ^  ^2^  ■, 

(7-27) 

which  implies 

(7-28) 

and 

9p{r,l3)  = 

2t(1  -  0t)  2(1  -  /5t)2  _  -2/5r  +  2/5V2  -  2(1  -  /5t)2 

(7-29) 

CO 

1 

CO 

1 

CM 

1 

Additionally, 

ifrh  -  fhr)  =  -2(4  +  4/5^2  -  /5(1  +  2r)2) 

(7-30) 

and 

iUh  -  fh^)  =  2r(l  -  2r)^ 

(7-31) 

are  needed.  Substituting  Eq.  (7-28),  Eq.  (7-29),  Eq.  (7-30),  and  Eq.  (7-31)  into 

Eq.  (7-4), 

gives 

4(-l  +  ;5T)2(-4-f /5  +  4/5r2)  ^ 

(7-32) 

Therefore, 

4 

l-f-4r2 

(7-33) 

and 

4t 

ol  = 

l-h4r2- 

(7-34) 

This  gives  a  relationship  beteen  a  and  j3^ 

i3  =  2-2yjl-  a2. 

(7-35) 
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7.3  TRANSIENT-POSITION-NOISE-REDUCTION  COST  FUNCTION 


To  determine  which  different  coefficient  relationship  are  achieved,  first  consider  the  transient- 
position-noise-reduction  cost  function,  Jntp'i  note  that 


/(T,/3)  =  2;3r2  +  (2-3r/3), 

(7-36) 

which  implies 

MT,/3)=4^r-3/3, 

(7-37) 

and 

//j(t,/3)  =  2t^  -  3r. 

(7-38) 

Also, 

(2  - /3r)(l  - /5t)2 

(7-39) 

which  implies 

9T{r,l3)  = 

2(-4  -1-  /3(1  t(4  -f  +  P{-2  +  t(-4  +  P  +  /3r)))))) 

(7-40) 

t2/32(4  -  2/3t  - /3)2 

and 

,  2{-l  +  Pr){8  +  p{-S  +  2ri-4:  +  p  +  l3r))) 

tI3^  (4  -  2/3t  -  /?)2 

(7-41) 

Additionally, 

ifrh  -  fhr)  =  -2(4  -F  4;0V2  -  /3(1  +  2t)2) 

(7-42) 

and 

(M-/M  =  2r(l-2r)2 

(7-43) 

are  needed. 

By  substituting  Eq.  (7-40),  Eq.  (7-41),  Eq.  (7-42),  and  Eq,  (7-43)  into  Eq.  (7-4), 

gives 

A{-l  +  pTf{-2  +  l3T{l  +  T))  ^ 

(7-44) 

Therefore, 

^  ~  t{1  +  t) 

(7-45) 

which  is  the  Benedict-Bordner  relationship. 
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7.4  TRANSIENT- VELOCITY-NOIS&REDUCTION  COST  FUNCTION 


To  determine  which  different  coefficient  relationship  are  achieved,  first  consider  the  transient— 
velocity-noise-reduction  cost  function,  Jnvt'-,  note  that 


/('T,^)  =J0, 

(7-46) 

which  implies 

/r(r,/3)  =  0, 

(7-47) 

and 

//3(t,/?)  =  1. 

(7-48) 

Also, 

-2-2/3(-l  +  T)r  +  A^ 

/3r(-4  +  /?  +  2/?r)  ’ 

(7-49) 

which  implies 

Pr(r,/3)  = 

-2/3(-l  +  r)  -  2/3r  -)-  2(-2  -  2/3(-l  +  r)r  -1-  0^t^) 

(7-50) 

/3r(-4-|-/3-|-2^T)  t(-4 -1- /3 -1- 2/3t)2 

-2-2/3(-l-t-r)T-|-/3V 

0t^{—4:  +  /3  +  2^t) 

and 

= 

-2(-l  +  T)r  +  2I3t^  (1  -1-  2r)(-2  -  2I3{-1  -1-  t)t  + 

(7-51) 

/5r(-4  +  ^  +  2Pt)  I3r{-A  +  l3  +  2/3r)2 

-2-2/3(-l-Hr)+/3V 

l3‘^r{-4  +  /3  +  2l3T)  ' 

Additionally, 

{frh  -  fhr)  =  /?(-4  -f  4tP  -t-  P) 

(7-52) 

and 

(fish  -  fh/s)  =  4t 

(7-53) 

are  needed.  By  substituting  Eq.  (7-50),  Eq.  (7-51),  Eq.  (7-52),  and  Eq.  (7-53)  into 

Eq.  (7-4), 

gives 

4-2Pt{1  +  t)) 

PT 

(7-54) 

Therefore, 

P=  ^ 

^  r(l-br) 

(7-55) 

which  is  the  Benedict-Bordner  relationship.  Thus,  we  have  the  interesting  observation  that 
either  transient  cost  function  leads  to  the  Benedict-Bordner  relationship.  It  is  an  interesting 
intellectual  challenge  to  determine  what  cost  function  leads  to  either  the  Kalata  or  to  the 
CTWN  relationship. 
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8  CONCLUSION 


The  performance  of  an  a  —  and  also  axi  cx  —  /3  —  ^  filter  for  several  deterministic  motion 
models  has  been  analyzed.  Closed  form  solutions  have  been  obtained  for  the  filter  performance 
statistics.  The  noise  reduction  ratio  is  used  as  a  new  means  of  determining  filter  lag.  Cost 
functionals  were  used  as  alternate  methods  for  determining  the  tracking  index.  This  report 
presents  a  general  solution  to  the  constant  gain  tracking  filters,  which  are  nothing  more  than 
matrix  difference  equations.  The  general  solution  can  be  used  to  compute  the  covariance  matrix 
imder  very  general  assiunptions  about  the  noise. 
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APPENDIX  A 
EVALUATING  SUMS 
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In  this  appendix,  we  evaluate  the  sums  that  occur  in  the  main  body  of  this  report.  The 
sums  are  of  the  form  of  sums  of  powers  of  sines  or  cosines.  We  use  the  complex  representation 
of  the  trigonometric  form  in  order  to  simplify  these  calculations. 

Let  y  =  re^®,  then 

(A-i) 

n=0  ^ 

if  -1  <  Re(y)  <  1.  This  can  be  used  to  evaluate  trigonometric  sums  by  noting 

OO 

lm{F{y))  =  S{r,9)  =  ^  sin(fc^)  (A-2) 

fc=0 


and 

OO 

Re{F{y))  =  C{r,0)  =  cos(/!:0) 

The  individual  sums  can  be  evaluated  by  applying  the  definition 


and 


/  SN  /  1  \  rsin(^) 

S(r.«)  =  Im(Fto))  =  Im  (3^)  =  +  = 

,  /  1  \  1— rcos(0)  1 

C(r,«)  =  E*(F(!/))  =  Re  ecs{9)  +  “  ‘ 


r  sin(^) 

~T~ 


—  r  cos{6) 
“0 


where 

/?  =  1  -  2r  cos(^)  +  r^. 


(A-3) 


(A-4) 


(A-5) 


(A-6) 


For  finite  sums,  define  the  sum  as 


F{y^k) 


k 


Ev' 


k 


=  i: 

n=0 


yk+l  _  1 

y-1 


(A-7) 


Hence,  we  have 

S(r,^)=Im(^(y))=Im 


T-fc+2  sin(fc^)  +  r  sin(^)  —  ^n{k  +  1)^ 

.  y-1  j 


(A-8) 
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and 


C(r,^)  =  Re(F(y))  =  Re 


''yfc+i  _  i\  cos{k)9  +  cos(fc  +  1)0  +  r  cos(0)  — 


y- 1 


The  first  sum  can  be  written  as  follows: 


(A-9) 


fc-1 


=  XI  sin(m0)  =  Im  [F{y,  fc)] , 

m=0 


(A-IO) 


=  Im 


iy‘‘  - 1) 
(v-i) 


=  Im 


rsin(0)  sin((fc  —  1)0)  —  r^sin(fc0) 

1  —  2rcos(0)  +  r2 


fc-l 


C(r,  0)  =  X  cos(m0)  =  Re  [F{y,  A:)] , 

m=0 


(A-ll) 


=  R,6 


(/-i) 


=  Re 


1  —  rcos(0)  +  cos{{k  —  1)0)  —  cos(A:0) 
1  —  2r  cos(0)  + 
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We  define  Si{r,9)  as 

fc-i 

Si(r,6)  =  ^2  'mr‘^sm{m6), 

m=0 


which  equals 


d 


k-l 


=  r—  r”^sin(m0) 
dr  ^ 


m=0 


Q-r  sin(g)  +  (g  +  /5fc)r*=+^  sin((fc  -  l)6i)  -  (a  -  (5{l  -  k))r‘ 
-  J2 


Similarly,  we  define  Ci  (r,  6)  as 


k-l 

Ci(r,9)  =  ^  mr^cos{m9) 

m—0 


which  equals 


o  fc— 1 

=  r—  r^cos{m9) 

-2r2  +  (2  -  a)r  cos{9)  +  [a  +  pk\  cos((fc  -  1)^)  -{a-  P{1 


(A-12) 

(A-13) 

sm{k9) 


(A-14) 

(A-15) 

—  A:)]r*cos(fc0) 
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We  define  S2  (r,  6)  as 


which  can  be  evaluated  as 


fc-i 


S2{t-,0)  —  ^  rr?r'^ sm{m6) 
m=0 


(A-16) 


(A-17) 


to  give 


=  {(2q:^  +  a/?  —  2j0)rsin(0) 


(A-18) 


+(2a^  —  2P  +  a/3  +  2a(3k  +  sin((A:  -  1)0) 


-(2a^  -  2/3  -  a/3  +  2af3k  +  /3^  -  2fc/3"  +  k^^y^  sin(fc0)}}\/3' 
Similarly,  we  define  C2(r,0)  as 

k-l 

C'2(r,0)  =  ^  mV’^cos(m0) 


which  can  be  evaluated  as 


to  give 


(A-19) 


(A-20) 


=  {— 4ar^  +  (4a  —  2a^  —  a/3)r  cos(0) 


(A-21) 


+(2a^  —  2/3  +  a/S  +  2a/3fe  +  cos{{k  —  1)^) 


— (2a^  —  2/3  —  a/3  4-  2a/3fc  +  0^  —  2k0^  +  k^0‘')r^  cos{k9)}\P' 
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